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We propose a new lattice framework to extract the relevant gluonic energy scale of QCD phenom- 
ena which is based on a "cut" on link variables in momentum space. This framework is expected to 
be broadly applicable to all lattice QCD calculations. Using this framework, we quantitatively de- 
termine the relevant energy scale of color confinement, through the analyses of the quark-antiquark 
potential and meson masses. The relevant energy scale of color confinement is found to be below 
1.5 GeV in the Landau gauge. In fact, the string tension is almost unchanged even after cutting 
off the high-momentum gluon component above 1.5 GeV. When the relevant low-energy region is 
cut, the quark-antiquark potential is approximately reduced to a Coulomb-like potential, and each 
meson becomes a quasi-free quark pair. As an analytical model calculation, we also investigate the 
dependence of the Richardson potential on the cut, and find the consistent behavior with the lattice 
result. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.Gc, 14.40.-n 
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I. INTRODUCTION 

In nature, there exist a large number of physical phe- 
nomena, and their scales range from the smallest scale 
(the Planck scale) to the largest scale (universe size). 
These phenomena are described in various ways depend- 
ing on their scales. For example, small-scale phenom- 
ena are described by elementary particle physics, nuclear 
physics, and so on, while large-scale phenomena by cos- 
mology, astrophysics, and so on. The scale characterizes 
each physical phenomenon and creates their hierarchy. In 
modern physics, the scale is one of the most fundamental 
and important concepts. 

In quantum chromodynamics (QCD) and some other 
field theories, there exists a characteristic phenomenon 
that describes the nontrivial appearance of a scale. In 
classical QCD Lagrangian, only dimensional quantities 
are quark masses, and therefore QCD is scale invariant 
in the chiral limit at the classical level. However, after the 
quantization, scale invariance is violated and an energy 
scale appears in a nontrivial manner. This is called as 
dimensional transmutation, through which many dimen- 
sional quantities are created in real QCD, such as "mass 
gap" , string tension, various condensates, and mass and 
size of hadrons. Dimensional transmutation is an impor- 
tant phenomenon in modern field theory P, Q ■ One of 
the most well-known energy scales in QCD is the QCD 
scale parameter Aqcd, which is defined through the run- 
ning coupling constant as{Q^) in perturbative QCD. It is 
the energy scale where the naive perturbative calculation 
does not work at all. 

The QCD running coupling constant depends on the 
energy scale as 
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where [}q = {llNc-2N f) / A&i:'^ , at one-loop level. There- 
fore, the behavior of QCD phenomena at high energy is 
entirely different from that at low energy. At high en- 



ergy or short distance, perturbative QCD is valid due to 
the asymptotic freedom [3, 'S]. In contrast, at low energy 
or long distance, the value of the QCD running coupling 
constant is large and nonperturbative effects are impor- 
tant. There exist many characteristic phenomena in non- 
perturbative region, such as quark confinement and chiral 
symmetry breaking d, d, 0, H, Hi • It is found that the 
appearance of the energy scale enriches the QCD phe- 
nomena. 

One important question is what energy scale is relevant 
for each QCD phenomenon. One may naively believe 
that the typical energy scale of nonperturbative QCD is 
Aqcd- However, Aqcd is just one scale that is defined 
through the running coupling constant. We do not know 
which energy components of gluons and quarks mainly 
contribute to a QCD phenomenon. Also, it is not trivial 
whether such relevant energy components are the same 
or different between different nonperturbative phenom- 
ena, for example, quark confinement and chiral symme- 
try breaking. 

In order to determine the relevant energy scale of QCD 
phenomena, we employ lattice QCD, which is the nonper- 
turbative and first-principle calculation of QCD [l^, [ll[ . 
In this paper, we propose a new framework in lattice 
QCD to determine the relevant gluonic energy scale sys - 
tematically, quantitatively, and nonperturbatively 
We study the relevant gluonic energy scale of color con- 
finement by this framework. 

The outline of the paper is as follows. In Sec. II, we 
formulate a lattice framework to determine the relevant 
energy scale for QCD phenomena. This framework is 
fairly general and readily applicable to actual lattice cal- 
culations. To investigate the relevant energy scale of color 
confinement, we calculate two fundamental examples of 
quenched lattice calculations; the quark-antiquark poten- 
tial in Sec. HI and meson masses in Sec. IV. In Sec. V, we 
perform an analytical study for the relevant energy scale 
of color confinement. We investigate the dependence of 
the confinement in the Richardson potential on the value 
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of the infrared cut, and compare it with the lattice re- 
sults presented in Sec. III. Finally, Sec. VI is devoted to 
summary and discussion. 



II. LATTICE FORMULATION 

The main concept of our framework is simple and gen- 
eral. It is to calculate a physical quantity, after artifi- 
cially cutting a certain region of momentum space. From 
the relationship between the cut region and the resulting 
quantity, we can determine the relevant energy scale of 
the quantity. This concept would be applicable to many 
theories in physics |13l] . Here, we introduce this concept 
to lattice QCD. The fundamental degree of freedom in 
lattice QCD is a link variable, which represents a gluon 
field on lattice. We consider a "cut" on link variables in 
momentum space, and investigate its effect on the result- 
ing quantity. 

The procedure for each gauge configuration is formu- 
lated as the following five steps. 

Step 1. We generate a gauge configuration by Monte 
Carlo simulation of lattice QCD under space-time peri- 
odic boundary conditions, and obtain a finite number of 
coordinate-space link variables 
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where a is the coordinate-space lattice spacing. Since the 
procedure is not gauge invariant, we fix the link variables 
with some gauge. It is desirable to choose a physically 
meaningful gauge on the lattice. In this paper, we mainly 
use the Landau gauge. The Landau gauge is well-known 
and frequently used also in continuum theory, and it gives 
a transparent connection between the link variable and 
the gauge field. 

Step 2. By carrying out discrete Fourier transforma- 
tion in the four-dimensional Euclid space, we define the 
momentum-space link variable U^{p) as 



(3) 



where TVsitc is the total number of lattice sites. The mo- 
mentum space represents a lattice with A'sito lattice sites, 
and boundary conditions are periodic. The first Brillouin 
zone of the momentum space is a four-dimensional hy- 
percube with each side of (— 7r/a, 7r/a]. The momentum- 
space lattice spacing is given by 
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where L is the number of lattice sites in each direction. 
The momentum-space lattice spacing corresponds to the 
minimum unit of momentum, and it has mass dimension. 

Step 3. We introduce a "cut" on Ufj.{p) in a certain 
region of momentum space. Outside the cut, U^{p) is 




FIG. 1: Two-dimensional schematic figure of the UV cut Auv 
and the IR cut Air on momentum-space lattice. The cut 
(shaded) region is described by Eq. The momentum- 

space link variable U^ijp) is replaced with the free variable 
U^il^'^ ij>) in the cut regions, ttp is the momentum-space lattice 
spacing. 



replaced by the free-field link variable 



Tip) = 1^ E 1 exp(iE.P.^.) = <^po. (5) 



U. 



Then, the momentum-space link variable with the cut is 
defined as 



U^{p) (inside cut) 

U^ip) ~ Spo (outside cut). 



(6) 



The concrete form of the cut can be taken arbitrarily. 
The most natural choice to estimate the energy scale is 



the cut by a four-momentum length 



jPf^Pt^^ 

which corresponds to a simple momentum cut in contin- 
uum theory. For example, an ultraviolet (UV) cut or an 
infrared (IR) cut by a four-momentum length is 



^ > Auv or < A 



IR- 



(7) 



These cuts are schematically depicted in Fig. [T] In this 
study, we mainly use these UV and IR cuts for numeri- 
cal calculations. Another possible choice is, for example, 
max(pi,p2,J53,P4) instead of which respects the lat- 
tice structure. Of course, we can cut not only the UV 
or IR region, but also intermediate-momentum region, 
anisotropic region, and so on. 

Step 4. To return to coordinate space, we carry out 
the inverse Fourier transformation as 



Ki^) - E ^l^(P) exp(-zE.P.a;,). 



(8) 
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Since this U'^{x) is not an SU(3) matrix, we project it 
onto an SU(3) element U^{x) by maximizing 

ReTV[t/^^x)tC/;(x)]. (9) 

Such a projection is often used in lattice QCD algorithms. 
By this projection, we obtain the coordinate-space link 
variable U^{x) with the cut, which is an SU(3) matrix 
and has the maximal overlap to U'^{x). 

Step 5. Using these link variables U^{x)^ we compute 
the expectation value of physical quantities in the exact 
same way as original lattice QCD. 

By repeating Steps 3-5 with different values of the cut, 
we can quantitatively determine the relevant energy scale 
of a physical quantity. Since we only have to replace 
U^{x) instead of Uf^{x), this framework can be applied 
to all lattice calculations. In addition, since the fast 
Fourier transformation is numerically easy task compared 
to gauge configuration generation or solver calculation, 
the computing time required for each cut is almost the 
same as that required for original lattice calculations. We 
expect that this framework can be broadly used to the 
analysis for the relevant energy scale of QCD phenomena. 

As for the gauge fixing, we mainly use the Landau 
gauge. In lattice QCD with the Euclidean metric, the 
Landau gauge is defined by the condition that globally 
maximizes the quantity 

F[[/] = ^^ReTr[/^(x), (10) 

X /I 

by the SU(3) gauge transformation. In terms of the gauge 
field Afj_(x), this condition is equivalent to minimizing 

J d'^xTilA^ix)^} (11) 

in Euclidean QCD, and it is a sufficient condition for the 
local condition df_,Afj,{x) — 0. Then, the gauge fluctua- 
tion is maximally suppressed, so that the link variable 
can be expanded as 

Uf,{x) = l + igaAf,{x) + --- (12) 

in a well-defined manner, and there exists a transpar- 
ent connection between the link variable and the gauge 
field. If one is interested in the direct information about 
the gauge field, it is possible to construct the similar 
framework with gAf^{x), which is renormalization-group 
invariant, instead of U^{x). For example, in the Landau 
gauge, gA^{x) is approximately extracted as 

gA^{x) = gA^{x) - |-TrA^(x), (13) 

9M^) ^ ^^{U,{x)-Ul{x)). (14) 

The trace part is subtracted for the traceless property of 
the gluon field. 



One may feel that our framework is similar to renor- 
malization. This framework is not equivalent to renor- 
malization transformation. In the case of renormaliza- 
tion transformation, the gauge field outside cutoff is in- 
tegrated out, however, in the case of our framework, the 
gauge field outside the cut is simply removed. Such a sim- 
ple removal of the gauge-field components breaks gauge 
invariance. Since the gauge transformation property is 
nonlocal and nontrivial in momentum space, gauge in- 
variance is lost when a certain momentum region is cut. 
If we do not fix the gauge, the resulting expectation value 
is zero, even in the case of a gauge-invariant operator. 
Then, gauge fixing is needed in the framework. Note, 
however, that the procedure itself is not restricted to a 
specific gauge, and we can use other gauges instead of the 
Landau gauge. By comparing the results obtained in sev- 
eral different gauges, we can check the gauge dependence 
of the relevant energy scale. 

In the case of a general gauge, since gauge fixing and 
momentum cut do not necessarily commute, it is nontriv- 
ial that these two conditions are satisfied simultaneously. 
In other words, U^{x) can deviate from the gauge-fixing 
condition which is originally imposed on [/^(x). However, 
in the Landau gauge, we have numerically checked that 
U^{x) almost satisfies the Landau gauge fixing condition, 
i.e., F[U^\ is maximized, even after the back projection 
to SU(3). 

III. INTERQUARK POTENTIAL 

We apply our framework to the analysis of the in- 
terquark potential in terms of gluonic energy compo- 
nents. The interquark potential is one of the most fun- 
damental quantities in QCD, and many types of in- 
terquark potentials have been studied by lattice QCD 
[3 [H, [M El- Among them, we here calculate the 
most basic potential, i.e., the quark-antiquark {QQ) po- 
tential. The QQ potential is extracted from the expecta- 
tion value of the Wilson loop, which is a gauge-invariant 
path-ordered product of link variables along a loop. 

A. QQ potential with UV/IR cut 

The QQ potential is expressed as a sum of one-gluon- 
exchange Coulomb potential and linear confinement po- 
tential as 

V{R)^(jR-^ + C, (15) 

where R is the distance between quark and antiquark. 
The physical value of the string tension a is approxi- 
mately 0.89 GeV/fm, and the Coulomb coefficient A is 
approximately 0.26. The constant C is physically irrele- 
vant, and its value depends on regularization. The QQ 
potential includes both perturbative and nonperturba- 
tive ingredients. In short range, it is dominated by the 
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TABLE I; Asymptotic string tension Uasym of the interquark potential with the IR cut Air. The resuhs obtained under four 
different conditions are hsted; the main result (with the four-dimensional hyperspherical cut in the Landau gauge) obtained in 
Sec. III-A, the result in the Coulomb gauge, the result with the hypercubic cut, and the result with three-dimensional formalism 
obtained in Sec. flf-C. The lattice QCD calculations are performed on 16* lattice with (3 = 6.0, where the lattice spacings are 
a ~ 0.10 fm and ap ~ 0.77 GeV. In original lattice QCD, the string tension aa^ is about 0.051. 



Am/ttp 


2 


Casymffl^ (Coulomb) 


o-asymfl^ (hypcrcubc) 


o"asyma^ (3-dim.) 


1.0 


0.0469(58) 


0.0289(58) 


0.0469(58) 


0.0433(51) 


1.1 


0.0311(49) 


0.0190(59) 




0.0198(43) 


1.5 


-0.0019(20) 


0.0024(25) 




-0.0034(9) 


2.0 


-0.0132(6) 


-0.0041(10) 


-0.0142(8) 


-0.0092(5) 


3.0 


0.0003(12) 


0.0058(10) 


0.0065(20) 


0.0024(13) 



perturbative one-gluon-exchange Coulomb potential. In 
long range, it is dominated by the linear confinement po- 
tential, which is purely a nonperturbative phenomenon. 
We determine the relevant gluonic energy scales of the 
Coulomb and confinement potentials. 

The numerical simulation in this subsection is per- 
formed with the isotropic plaquette gauge action with 
(3 = 6.0. The lattice size is 16^ and periodic bound- 
ary condition is imposed, and then the correspond- 
ing momentum-space lattice is also a 16^ isotropic lat- 
tice. The gauge configuration number is 50, however, 
for statistical improvement, we average all the parallel- 
translated Wilson loops in one configuration and apply 
the APE smearing method [l3|. The coordinate-space 
lattice spacing a is about 0.10 fm, and the lattice vol- 
ume is (1.6 fm)"*. The momentum-space lattice spacing 
ap is about 0.77 GeV, and the momentum-space lattice 
volume is (12 GeV)^. In this paper, a dimensional quan- 
tity in coordinate space is scaled with a by the standard 
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FIG. 2: The QQ potential with the UV cut Auv plotted 
against the interquark distance R. The lattice QCD calcu- 
lation is performed on 16* lattice with (3 — 6.0. The unit is 
scaled with the lattice spacing, a ~ 0.10 fm or Up ~ 0.77 GeV. 
The broken line is the original QQ potential in lattice QCD. 



lattice convention. Similarly, a dimensional quantity in 
momentum space, such as the cut Auv or Air, is scaled 
with ttp. 

First, we introduce the UV cut Auv to the QQ poten- 
tial. The resulting potential is shown in Fig. [2l Since 
the maximum value of the four-momentum length is 
\fv^ jap = 16 in our lattice, the standard lattice result 
(broken line) corresponds to the case of Ayv/op = 16- 
When the value of the UV cut becomes smaller, the 
short-range Coulomb potential gradually decreases. The 
constant term C also decreases by the UV cut. This is 
because the constant term is mainly given by the lat- 
tice regularization for the UV singularity of the Coulomb 
potential. On the other hand, the long-range linear po- 
tential is almost unaffected by the UV cut. When we 
fit the QQ potential by Eq. (llSp . the string tension cr is 
insensitive to the value of the UV cut. At Auv/op = 2, 
the QQ potential becomes only the linear potential, and 
its fitting result is aa^ = 0.0459(15), A = -0.056(5), 
and Ca — —0.095(7). The string tension differs slightly 
from its original value aa^ ~ 0.051, and the Coulomb 
coefficient and the constant term are nearly zero. 
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FIG. 3: The QQ potential V{R) with the IR cut Air. 
notation is the same as Fig. (2] 



The 
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Next, we show the QQ potential with the IR cut Air in 
Fig. [31 and Hst the value of the asymptotic string tension 
Casym in the secoud column of Table HI The asymptotic 
string tension (Jasym is estimated by fitting the QQ poten- 
tial in 3 < R/a < 9 with a linear function tTasymfi+ const. 
In contrast to the UV case, the long-range linear poten- 
tial is affected by the IR cut. At Aift/ap = 1 and 1.1, the 
string tension decreases slightly from its original value 
a ~ 0.89 GeV/fm. In Aia/ap > 1.5, the long-range lin- 
ear potential disappears, and the QQ potential becomes 
the perturbative Coulomb potential. Because Op ~ 0.77 
GeV, the physical energy scale for the vanishing of color 
confinement is about 1 GeV. From Fig.[3]and TableHl we 
observe that the asymptotic string tension is a small neg- 
ative value at Am/up — 1.5 and 2. This would suggest 
that, although the confinement potential asymptotically 
disappears, it slightly survives only in the intermediate 
range. At AiR/op = 3, the QQ potential becomes com- 
pletely flat in R/a > 3. 

From these results, we conclude as the following. The 
perturbative and nonperturbative parts of the QQ poten- 
tial are decoupled in the momentum space of the gluon. 
The relevant energy scale of confinement is below about 
1 GeV in the Landau gauge. If the gluon is restricted 
to this energy region, the QQ potential becomes only 
a linear potential. In contrast, if this gluon is cut, a 
Coulomb-like potential is obtained. 



is restricted to discrete 



/y^P^PM = 0, ap,\/2ap,\/3ap,- 



(16) 



Then, we can only take discrete variation on the value of 
the cut. For example, the lattice calculations in the range 
< Air < Op yield the same result. This is a kind of dis- 
cretization error in momentum space. In order to achieve 
a finer resolution in momentum, we must calculate with 
a smaller momentum-space lattice spacing, i.e., a larger 



TABLE II: Lattice spacings with different lattice couplings 
P = 2Nc/g'^. The coordinate-space lattice spacing a and 
momentum-space lattice spacing Up are listed. 



/3 


lattice size 


a [fm] 


Up [GeV] 


5.7 


16* 


0.19 


0.41 


5.8 


16* 


0.14 


0.55 


6.0 


16* 


0.10 


0.77 
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FIG. 4: The Auv-dependence of the asymptotic string tension 
o"asym iu physical unit. The vertical error bar is the standard 
statistical error, and the horizontal error bar is the range that 
yields the same result due to the discrete momentum. The 
original value of the string tension is cr ~ 0.89 GeV/fm (bro- 
ken line). 
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B. More quantitative estimate for the energy scale ^ 
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For a more quantitative argument, we need higher ac- 
curacy in momentum space. Since the minimum momen- 
tum on lattice is the momentum-space lattice spacing Op, 
the four-momentum length 
values as 



p=6.0 1— e— 1 ■ 

P-=5.8-i L 
p=5.7 1— 1 



0.5 



1 1.5 
A|R [GeV] 



FIG 
in Fig. H 



5: The AiR-dependence of the asymptotic string tension 
asym iu physical unit. The notation is the same as that used 



coordinate-space lattice volume. In addition, it is desir- 
able to calculate with several different momentum-space 
lattice spacings. 

For more quantitative estimate on the relevant energy 
scale of color confinement, we calculate the QQ poten- 
tial on 16^ lattice with /? — 5.7, 5.8, and 6.0. As listed 
in Table |lTl the corresponding momentum-space lattice 
spacings are Up ~ 0.41 GeV, 0.55 GeV, and 0.77 GeV, 
respectively. The coordinate-space lattice spacings are 
determined so as to reproduce the string tension a to be 
0.89 GeV/fm. Other conditions are the same as before. 

In Fig. m we show the Auv-dependence of the asymp- 
totic string tension (Jasym • The asymptotic string tension 
Casym IS estimated by fitting the QQ potential with a lin- 
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tion for the Coulomb gauge is to maximize the quantity 
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FIG. 6: Another example of a cut on momentum-space lat- 
tice. The cut (shaded) region is described by Eq. (|18|l . Other 
notations are the same as those used in Fig. [T] 



ear function CTasym^ + const, in 0.3 fm < i? < 0.9 fm. In 
this figure, while the error bar in the vertical direction 
represents the standard statistical error, the error bar in 
the horizontal direction represents not the statistical er- 
ror but the range which yields the same result due to 
the discrete momentum. The asymptotic string tension 
is almost unaffected in Auv > 1-5 GeV, and its value 
is unchanged from its original value a ~ 0.89 GeV/fm. 
In Auv < 1-5 GeV, the value of the asymptotic string 
tension significantly decreases. 

We analyze the IR case in the same way and show the 
result in Fig. [5l All the results consistently suggest that 
the energy scale of the vanishing of the asymptotic string 
tension is around Air = 1.2 GeV. Above this energy 
scale, the asymptotic string tension is almost zero and 
the confinement potential asymptotically disappears . 

From these quantitative analyses, we conclude that 
the relevant energy scale of color confinement is below 
1.5 GeV, i.e., color confinement originates from the low- 
energy gluon components below 1.5 GeV. 



C. Consistency check with other conditions 



^^ReTrC/,(x). 

X j = l 



(17) 



The values of iTasym with the Coulomb gauge are listed 
in the third column of Table U The value itself depends 
on the gauge choice, however, the energy scale for the 
vanishing of confinement is the same, that is, about 1 
GeV. 

(ii) Our framework has an ambiguity of the cut form. 
For example, since the rotational invariance is broken on 
lattice, it is worth introducing a cut which respects the 
lattice structure of momentum space. Then we calculate 
with the cut by max(pi,p2,P3,P4) instead of y^p^, i.e., 

max(pi,p2,P3,P4) > Auv or max(pi,p2,P3,P4) < Air, 

.(18) 

instead of Eq. This cut forms a four-dimensional 

hypercube, as shown in Fig. [51 The values of (Tasym with 
this cut are listed in the fourth column of Table HI The 
result estimated by max(pi,p2,P3,P4) is consistent with 
that estimated by ^/p^ ■ 

(iii) Our framework can be easily extended to the 
spatial three-dimensional formalism. The spatial three- 
dimensional formalism is related to the Richardson po- 
tential calculation, which is shown in Sec. V. We only 
have to perform the three-dimensional Fourier transfor- 
mation and consider a cut by the three-momentum length 
\p\. The result of the three-dimensional formalism is 
listed in the fifth column of TableU In the Landau gauge, 
the relevant three-dimensional momentum scale of con- 
finement is found to be below about 1 GeV, which is 
almost the same as the relevant four-dimensional energy 
scale. 



IV. MESON MASSES 

In this section, we apply our framework to the analysis 
of meson masses in quenched lattice QCD. We calculate 
the mass of a pion (pseudo-sealer meson) and a p-meson 
(vector meson) using two different fermion actions, i.e., 
the clover fermion action and the staggered fermion ac- 
tion. 



We have calculated with other conditions for the check 
of consistency: (i) a different gauge in Step 1, (ii) a differ- 
ent cut in Step 3, and (iii) the three-dimensional Fourier 
transformation. The results obtained under these three 
conditions are consistent with the previous results. Here, 
we briefly summarize these calculations and results. 

(i) Since our framework is not gauge invariant, it is 
important to check that the resulting energy scale does 
not depend drastically on the gauge choice. We calculate 
with the Coulomb gauge for the gauge choice in the Step 
1, instead of the Landau gauge. The gauge- fixing condi- 



A. Clover fermion 

The clover fermion action is an 0(a)-improved Wilson 
fermion action, and the quark mass is described in terms 
of the mean-field-improved hopping parameter n ^^^^ 
[iO, HH. For the details of the improvement, see Ref. [l7l |. 
The Wilson fermion is the simplest fermion on lattice, 
but it explicitly breaks chiral symmetry for avoiding the 
doubling problem. 

The numerical simulation is performed on 16''' lattice 
with (3 = 6.0. The lattice spacings are the same as before, 
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FIG. 7: The Auv-dependence of pioii mass and p-meson 
mass nip. The quark propagator is calculated by the clover 
fermion action with the hopping parameter k. The data at 
Auv/flp = 16 is the standard lattice result. The unit is scaled 
with the lattice spacing, a ~ 0.10 fm or ap ~ 0.77 GeV. 



FIG. 9: The AiR-dependence of meson masses with the stag- 
gered fermion. The bare current quark mass is taken as 
ma = 0.01, 0.02, 0.06, and 0.10. The unit is scaled with 
the lattice spacing, a ~ 0.10 fm or ap ~ 0.77 GeV. 
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and the p-meson mass decrease by removing the IR gluon. 
In addition, in contrast to the case of the UV cut, these 
meson masses degenerate in Am/ap > 2. In the calcula- 
tion of the interquark potential, the confinement poten- 
tial vanishes in this energy region of Air. Therefore, this 
degeneracy in meson masses suggests that the quark and 
antiquark in such mesons become unbound or, if possible, 
very narrowly bound. This state is called as "quasi-free" , 
which implies that, in this state, the Coulomb interaction 
exists, but the confinement potential does not exist. 



FIG. 8: The AiR-dependence of meson masses with the clover 
fermion. The data at Am/ap = is the standard lattice 
result. The notation is the same as Fig. [T] 



i.e., a ~ 0.10 fm and ap ~ 0.77 GeV. Three different 
values of hopping parameters are taken as k = 0.1200, 
0.1300, and 0.1340, and the corresponding pion masses 
are about 2.9 GeV, 1.8 GeV, and 1.3 GeV, respectively. 
The configuration number is 100 here. 

Figure [7] shows the pion mass to^ and the p-meson 
mass TUp with the UV cut. In this figure, the data at the 
far right side, i.e., Auv/op = 16, is the standard lattice 
result. As the value of the UV cut becomes smaller, both 
the pion mass and the p-meson mass gradually decreases. 
This decrease corresponds to the mass generated by the 
UV gluon which dresses the quarks in the mesons. 

Figure [S] shows TOtt and with the IR cut. The data 
at the far left side, i.e., Auv/op = 0, is the standard lat- 
tice result. As in the case of the UV cut, the pion mass 



B. Staggered fermion 

The staggered fermion is another well-known fermion 
in lattice QCD, which has a kind of chiral symmetry. 
This fermion is often used for investigating the chiral 
property of hadrons [13, ■ 

The staggered fermion calculation is performed on 
IG'^ X 32 lattice with (3 = 6.0. Then, in the spatial direc- 
tion, the momentum-space lattice spacing is Up ~ 0.77 
GeV, and in the temporal direction, the momentum- 
space lattice spacing is ap/2. The bare current quark 
mass m in the staggered fermion action is taken as 
ma = 0.01, 0.02, 0.06. 0.10, and the corresponding pion 
masses are about 0.48 GeV, 0.67 GeV, 1.2 GeV, and 1.5 
GeV, respectively. 

The dependence of the meson masses on the IR cut 
is shown in Fig. [51 The staggered fermion exhibits the 
same behavior as the clover fermion, i.e., both m,r and rrip 
decrease by the IR cut, and they degenerate in Ain/op > 
2. 

Next, we consider the dependence of the meson masses 
on the quark mass. For chiral extrapolation, the squared 
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meson masses are fitted with quadratic functions as 

(m^a)^ or [rapof = C2{ma)'^ + Ci(ma) + Cg, (19) 

where C2, ci, and co are fitting parameters. In the stan- 
dard lattice calculation, the masses of p-meson and many 
other mesons are linear functions of the quark mass, and 
they are finite value even in the chiral limit m ^ 0. On 
the other hand, the pion exhibits characteristic proper- 
ties. It obeys the well-known relation near the chiral 
limit, 

flml^-m{qq), (20) 

which is called the Gell-Mann-Oaks-Renner relation. 
Here, {qq){— {uu + dd)) is the two-flavor quark conden- 
sate. In the chiral limit, the pion becomes massless be- 
cause it is the Nambu-Goldstone boson associated with 
chiral symmetry breaking. 

We show the quark mass dependence of the pion mass 
in Fig.[Tni and that of the p- meson mass in Fig. [11] When 
the IR cut is introduced, the p-meson mass uniformly 
decreases with the linear extrapolation form unchanged. 
In contrast, the pion loses its character in Am/ap > 2, 
where the pion mass degenerates to the p-meson mass. 
The extrapolation function becomes a linear function. In 
other words, the Gell-Mann-Oaks-Renner relation is bro- 
ken, and the pion is no longer the Nambu-Goldstone bo- 
son. We note that this behavior of the pion originates 
from the vanishing of confinement, rather than chiral 
symmetry restoration. In Am/ap > 2, the pion is not 
a bound state and far from the Nambu-Goldstone boson. 

Chiral symmetry breaking is another important topic 
in nonperturbative QCD. It is interesting to accurately 
investigate chiral properties with a small quark mass and 
unquenched calculation. 

V. ANALYTICAL MODEL CALCULATION 




0.05 0.1 



ma 

FIG. 10: The rho mass nip plotted against the bare current 
quark mass m in the staggered fermion action. Chiral extrap- 
olation is performed with {nipa)'^ — C2(ma)^ + ci{ma) + cq. 
The data with Air — (broken line) is the standard lattice 
result. The unit is scaled with the lattice spacing, a ~ 0.10 
fm or ap ~ 0.77 GeV. 
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A. Richardson potential 

From the lattice QCD results, it is found that the 
relevant energy scale of color confinement is below 1.5 
GeV. In this section, we introduce the IR cut to the phe- 
nomenological interquark potential which includes the 
confinement potential, and compare it with the lattice 
result. For this purpose, we analyze the Richardson po- 
tential dii. 

The Richardson potential is the phenomenological in- 
terquark potential constructed so as to reproduce the 
Coulomb plus linear structure. This potential is defined 
by the one-dressed-gluon-exchange amplitude which is 
proportional to 



V{p^) = -Ci 



where 



1 



/3oln(l+pVA2 



(21) 



(22) 



FIG. 11: The pion mass m-^ plotted against the bare current 
quark mass m in the staggered fermion action. The notation 
is the same as Fig. [TOl 



Here, Cp = 4/3 and (3o = (H^c - 27V/)/487r2. A is the 
only parameter in this model. This coupling is similar 
to the standard QCD coupling, except for 1 in the ar- 
gument of logarithm. After the integration of the time 
component, the Richardson potential is obtained by the 
three-dimensional Fourier transformation as 



V{R) = 



d^p 



(27r)3 



JR-P 



V{f) 



where 



87r/3o 



dt- 



A'R- 



/(Ai?) 



R 



R 



1 



t [ln(t2 - 1)]2 +7r2 



(23) 



(24) 



9 



We set Nf = to compare this potential with the QQ 
potential in quenched lattice QCD, and A = 0.48 GeV so 
that the coefficient of the linear potential is equal to the 
physical string tension a ~ 0.89 GeV/fm. 



B. Richardson potential with IR cut 

We introduce the IR cut Air to the Richardson poten- 
tial. We consider following two ways for the IR cut. 

One is the simple IR cut by the three-momentum 
length on the Fourier transformation, 



JR-p 



'b1>AiR (2^)^ 
The other is the change of the functional form as 



V{p 



7)2 -I- A 2 ■ 
P + ■'MR 



(25) 



(26) 



The advantage of the latter way is that we can analyt- 
ically calculate the momentum integral. The analytical 
derivation is given in Appendix. The result is 



V{R)^ 



where 



87r/3o 



-XAR 



fxjAR) 
R 



In(A^-l) 
[ln(A2-l)]2+7r2 



(A>1) 



RT^ (1>A>0) 



h{x) = AP 



dt 



te- 



1 



t2-A2 [In(i2-l)]2 + ;r2 



,(27) 

(28) 
■ (29) 



and A = Air /A. The symbol P indicates the principal 
value of the integral, which is necessary in the case of A > 
1. This functional form becomes the original Richardson 
potential in the limit of Air — > 0, apart from an irrelevant 
constant. 

These two ways to introduce the IR cut yield almost 
the same results. The result of the former way is shown in 
Fig. [T^ The irrelevant constant is arbitrarily subtracted 
in the figure. As in the case of the lattice result, the 
string tension decreases when the IR cut is introduced. 
In Air > 1 GeV, the linear confinement potential disap- 
pears and the interquark potential becomes a short-range 
Coulomb-like potential. This behavior is consistent with 
our lattice result. While the Richardson potential is only 
a phenomenological model and its confinement potential 
is set by hand, it well reproduces the AiR-dependence of 
the interquark potential obtained by lattice QCD. 



VI. SUMMARY AND DISCUSSION 

We have formulated the lattice framework to study the 
relevant energy scale of QCD phenomena. We cut link 



U'.:y A|F,=0 

V A|F|=oSGeV 

A|R=0.7GeV 

A|Fi=1.0GeV 

A|F|=1.3GeV 

0.2 0.4 0.6 0. 

R [fm] 



FIG. 12: The Richardson potential with the IR cut Air. The 
solid line is the original Richardson potential (|23|l . The irrel- 
evant constant is arbitrarily subtracted. 



variables in a certain region in momentum space, and cal- 
culate the expectation value of physical quantities with 
the link variables with the cut. Using this framework, we 
can determine the giuonic relevant energy scale of QCD 
quantities. This framework is broadly applicable for all 
lattice QCD calculations. 

From the asymptotic string tension of the QQ po- 
tential, we have found that the relevant giuonic energy 
scale of color confinement is below 1.5 GeV. In fact, 
the string tension is almost unchanged even after cut- 
ting off the high- momentum gluon component above 1.5 
GeV. When we cut the low-momentum component be- 
low 1.2 GeV, the confinement potential asymptotically 
disappears. Then, when we cut the low-momentum com- 
ponent up to 1.5 GeV, almost all the contributions to 
the string tension disappear in the whole range of the 
interquark distance. In addition, by cutting the infrared 
gluon in this energy region, the two quarks in mesons are 
not confined and become quasi-free. 

It is often believed that the typical energy scale of 
nonperturbative QCD is Aqcd- However, the value of 
the QCD running coupling constant is large even above 
Aqcd- In fact, our result suggests that the relevant en- 
ergy scale of color confinement is larger than Aqcd • This 
fact would be useful as a reference for developing effec- 
tive theories. For example, the relevant energy scale de- 
termines the cutoff in a low-energy effective model. Al- 
though the value of the cutoff should be based on some 
physical reasoning, its microscopic derivation is difficult 
in many cases. Since lattice QCD is the first-principle 
calculation in QCD, the relevant energy scale obtained 
here would provides a strong physical reasoning for the 
cutoff value. If we set the ultraviolet cutoff to be 1.5 GeV 
in momentum integral, we can safely pick up the contri- 
bution to the confinement. Further, the relevant energy 
scale determines degrees of freedom which appear in the 
effective field theory |25|]. 
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The relevant energy scales can differ even among many 
nonperturbative phenomena. Our framework is consid- 
ered to be a useful tool for understanding the energy 
scales of many QCD phenomena. 



For performing the integration in Eq. (jA2p , we consider a 
contour integral in the complex g-plane shown in Fig. 1131 
The radial part drops with an infinitely large radius. The 
two pole contributions at q = and iX are given by 
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APPENDIX A: ANALYTICAL DERIVATION 

In this section, we derive Eq. (j27p from the Fourier 
transformation of Eq. 



V{R) = - 



Cf f d^P 



iR-p 



(Al) 



/3o J (2^)3 (p2 + A2p) ln(l + p2/A2) 
After performing the angular integration, we obtain 

r-i A POO pp'ii^Q 

V{R) = Im / dq—^ -If— ^,(A2) 



2tt^Po r Jo {q^ + A^) ln(l + q^) ' 



where 



Ai?, q = |p1/A, a = AiR,/A. 



(A3) 



hx = 



ITT 1 

y A2 



2 ln(A2-l)+iTr -*-) 

TRT^ (1>A>0), 



(A4) 
(A5) 



respectively. The integral along the imaginary axis is 
given as 



Jo (g2 + A2)ln(l+g2) 

A * (x^ - t^)[Ht^ - 1) + 17T] 

-P I dt — — ^T^^, —. (A6) 







(A2-<2)in(l-t2)- 



The symbol P indicates the principal value of the inte- 
gral. By inserting these, we get 



V{R) 



Cf Ar 2 



87r/3o r 



— -2hxe-'^ + hir) , (AT) 



where h\ and f\{r) are given as Eq. ([28]) and Eq. ([29]) . 
respectively. 
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